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INTRODUCTION 


'  Knowledge  of  the  errors  in  a  diffraction  measurement  of  residual  strains 
and  atrassea  is  useful  information,  not  only  in  its  own  right,  but  also  be¬ 
cause  it  permits  automation  of  a  measurement  to  an  operator  specified  pre¬ 
cision^  ~There  are  three  sources  of  thase  errors: 

(  1)  Instrumental  effects:,  primarily  due  to  sample  displacement,  sep- 


ration  of  the  6  and  20  axes  of  the  diffractometer,  and  beam 

2 

divergence.  All  three  can  be  estimated  ,  or  minimized  by  em- 

3 

ploying  parallel  beam  geoawtry. 


(2) 


Qsv  d  [3) 


Uncertainties  in  x-ray  elastic  constants;  which  can  now  be  eval- 

4  — ' 

uated.  - - "'T" 

v  (t°  /*?*) 

Errors  ii  liffractlon  peak  position  related  to  counting  sta¬ 


tistics.  '"equations  to  evaluate  this  source  have  been  developed 


in  Ref.  1  for  the  case  of  a  stress  state  for  which  all  (1 
■  1,2,3)  ■  0,  with  the  direction  "3"  normal  to  the  sample  surface, 
see  Fig.  1.  This  means  that  the  stresses  lie  only  in  the  sur¬ 


face,  e.g.,  a  biaxial  stress  state  <^a 

<*a  <%s 


There  are,  however, 


numerous  situations  when  the  normal  components  are  appreciable 
in  an  x-ray  measurement^ ’ ^  and  this  is  generally  the  case  for 
neutron  diffraction  because  with  neutrons  the  beam  can  sample  a 
sizeable  volume,  at  a  significant  depth  below  the  surface^.  It 


FIG.  1:  The  axial  system.  Strains  are  measured  with  diffraction  by 
.  measuring  the  change  In  spacing  of  planes  normal  to  the  I* 
direction.  (The  axes  P^  define  the  specimen  surface.) 


TABLE  I:  STRESS  TENSORS  (  AND  STANDARD  DEVIATIONS) 
FOR  SPECIMEN  C3,  REF.  5 


DATA  SET  1 


537.62 

(161.94) 

-24.03 

(78.81) 

-39.15 

(4.58) 

550.04 

(161.66) 

2.31 

(3.56) 

78.29 

(130.57) 

DATA  SET  2 

520.60 

(137.25) 

-4.03 

(66.60) 

-34.17 

(3.21) 

555.19 

(137.03) 

0.11 

(2.69) 

82.20 

(110.67) 

DATA  SET  3 

535.03 

(158.28) 

-20.13 

(77.06) 

-40.19 

(5.72) 

555.98 

(157.99) 

-0.98 

(4.56) 

86.66 

(127.69) 

DATA  SET  4 

538.53 

(146.23) 

-30.63 

(70.95) 

-38.03 

(3.83) 

565.37 

(146.14) 

0.76 

(3.89) 

88.18 

(117.92) 

AVERAGE 

REFERENCE  5 

532.95 

-19.69 

-37.89 

541  -20 

-38 

556.65 

0.55 

565 

1 

83.83 

86 

values  given  In  MPa;  V(<^  )  ■  0.00016  k 


3 


is  the  purpose  of  this  paper  to  derive  equations  to  evaluate  the  counting 
statistical  error  for  the  entire  three  dimensional  strain  (or  stress)  tensor, 


<h  1 

<h> 

<**3 

<h* 

0*3 

BASIC  EQUATIONS 

He  begin  with  the  general  equation  for  the  strains  (e^ )  and  how  these 
affect  the  interplanar  spacing  "d".  (Refer  to  Fig.  1  for  the  axial  system.) 
The  measurement  is  made  in  the  A  direction,  with  a  sample  tilted  A  from  the 
normal  position  (which  is  with  the  surface  normal  bisecting  incident  and 
scattered  beams).  Primed  quantities  refer  to  strains  in  the  co-ordinate 
system,  unprimed  terms  are  in  the  P^  system. 

<«*',  ■  (d^  -  (]«  )/do  ■  C<*ii>  cos*#  +  <tt*>  sin?#  +  <elt>  sin2# 

-<«j s  >3  sin?#  +  <*j j  >  +  [<el,  >  cos#  +  <**,  >  sln#]sln2#  (1) 

Note  that  the  stress-free  spacing,  da  ,  is  involved.  While  this  term 
can  be  eliminated  for  a  biaxial  stress  state,  this  is  not  possible  for  a 
general  strain  or  stress  tensor,  and  the  reader  may  consult  Ref.  8  for  a 
discussion  of  problems  associated  with  the  measurement  of  this  quantity. 

When  q,  or  f  0,  is  not  linear  with  sin?#  and  has  different  curvature 
for  •»#  and  -#.  The  carats  Imply  that  the  strain  values  are  averaged  over  the 
depth  of  penetration  of  the  incident  x-ray  (neutron)  beam  and  this  is  to  be 
understood  in  what  follows,  as  this  additional  notation  is  eliminated  below. 

Next,  we  define  terms  which  Involve  measurements  of  d.  .  at  plus  and 

5  *»▼ 

minus  #  tilts  of  the  surface  normal. 

-  =  wa**.  *  *^]  -  <V  ♦  dH-)/2d>  w 

a  t„  +  [cxl  cos1#  ♦  Cat  ■it'?#  +  <!t  «in2#  -  Cja])sin?#  (2a) 

Clearly,  ^  ,  should  be  linear  with  sin?#  and  CjS  is  the  intercept,  regardless 
of  #. 


a,  =  -  d#  ^.)/2d, 

•  Ccxa  cos#  +  Cat  sln#]sln|  2#  | . 


(2b) 


Therefore,  %  is  linear  vs.  sin]24*| . 


'  4 


Let: 

h  -  d^  /dsin*4., 
it  -  deg  /dsin  1 2<|>  | 

Then,  et:  •  -  0* ,  a  ■  *1X  -  , 

•  ■  90*,  ■  Cat  “  *j*» 

•  "  45* ,  4s  +  **»  )+*ia  “  *»»  > 

-  Caa  +  1/2  (0  lx  +  to  ft  >• 


end  similarly: 

at  •  -  0*:  ,1,  -  exl  , 
at  *  ■  90* :  aa  fa  ■  tai  • 


(3a) 

(3b) 


(3c) 


(3d)  \ 


Knowledge  of  the  strain  tensor  permits  the  calculation  of  the  stress  com¬ 
ponents  (°jj)  from: 


otJ  -  [l/2Sa  (hki)]”1  {clj  -  (hkl)/[3%  (hkl) 

♦  1/2S,  (hk£)]}«  [exx  +  c*a  +  s,,]}. 


(4) 


Here,  6^  is  the  Kronecker  delta  function  and  and  1/ 2S,  are  the  x-ray 
elastic  constants  which  depend  on  the  indices  of  the  diffraction  peak,  hki.  . 
(For  an  isotropic  solid  these  values  are  -v/E  and  (1  ♦  v)/E  respectively.) 


VARIANCES  DDE  TO  COUNTING  STATISTICS 


For  a  function  X  ■  f(^  ,  x,  ,  x, ...),  assuming  the  x  are  independent,  the 

9  n 

variance  (V)  is: 


V(X)  -  (£|  )2  V(*  )  ♦  (|2  )2  V(*  )  ♦  (||  )2  V(*  )  + 


(5) 


For  the  straight  iine,  y  ■  mx  ♦  b,  the  slope  and  intercept  is  given  by: 

l  <*t  ■  *>  (yt  -  y) 


Z  (x.  -  x)* 
l  x 


b  -  (Zyt  -  mXx^/N 


(6a) 

(6b) 


where  N  Is  the  msaber  of  data  points. 


Employing  Eq.  (5): 


\  (yi  "  *  [f  (x<  *  * 

- - i -  V(x)  +  i - i - 

I  (x  -  iV  I  (x.  -  x 

1  i  1 


Therefore : 


V(b)  . 


ZCXj^  -  x)a 


Kx  -  iy*  Ky  - 
.  V(m) - J -  .  {  - i— 


■  ?(X1 

KXj^  -  x)a 


V(y)} 


Therefore;  in  terms  of  %  vs.  sirf^: 


Z  (sii^4»  -  siir^j 
i  1 


V(sirf*) 


,r^-w»T.. . 

- -  V(*  ) 

zCsirf  if  -  sirf^)* 
i  1 


The  variance  in  *  is  negligible,  so  the  first  term  can  be  ignored. 


Also,  from  Eq.  (2a): 


v(°-  wfcj*  fe 


v(d„.)  •*•  373^7  v<^)  (8) 


+  Writing  Bragg's  law  in  the  form  d  -  adopting  the  convention  that 

0  ,  0  are  the  0  values  (in  degrees)  for  the  peaks  at  +$,  -4,  respectively, 
and  employing  Eq.  (5): 

■  («/180)* (Acosf+/2sirf e+)1  V(20+)/2  (9) 

and  similarly  for  V(d#  ^_).  Recalling  Eq.  (2a): 


5^  '  -  zap-. 


(10a) 


~ ^  V  /  • 


a 


Thus,  we  may  rewrite  Eq.  (7): 

|(aiiP<|>1  -  sin*4») 


Vdi  ) 


Z(siiP4»^  -  iI?T)a| 


•jrr% — { («/180)*  ~  [(££S^_  )*v  (2e+) 
**•  8  sitPe+  1 


+  (£22®-  )a  v.  (2e“)]  +  (df* /<$)v (d„)} 
sin1  e“  1  + 


(ID 


In  a  similar  manner  for  a,  vs.  sin|2t|i|,  where  a,  =  (d^  -  d#^_)/2d°  «8~/2d°  : 

|sln|24>1|  -  ;din|24.|) 


Z(sin|2*  |  -  sinj24,|)2 

a  1 


2  ./ 


+  (£2Si  )»  V .  ( 28“ ) ]  +  (d*/dg)  V(d,  )} 

sin*  0~ 


We  now  propagate  these  values  into  the  strain  and  stress  tensors. 
THE  STRAIN  TENSOR 

Abbreviating  the  Intercept  of  4  vs.  sin* 4.  as  I,  then  at  any  t: 

eS3  -  I  (of  %  vs.  sin*4), 

V(e,  j  )  -  V(ix  )  +  V(I) 


(12) 


(13a) 

(13b) 


V(I) 


Z^iiP^  -  sirP^)*  ZtsirP^  -  sin2 <|i) 
N  Z(sin?4»i  -  sirP^) 


V(*  i> 


.  [z(aiiP4»  -  sirP  <Ji)? 

—  - i -  V(s1 

N  Z(sin*4»i  -  airP^)* 


i) 


(13c) 


Now,  from  Eqs.  (3c),  at  *  -  0* : 


#  fi  ■  *11  *  *u  "*n  "I*  (14) 

V(elx  )  -  2V(,  fx  )  +  V(I). 

Similarly,  for  •  -  90* : 


V(*aa)  -  VC,,!*  )  +  V(,  ^  )  ♦  V(I),  (15) 

and  for  •  -  45* : 

*i»  ■  **  fi  +  *»»  “  0*5  (ei4  +  «»a  ) 

+  *4  i  -  0.5  (,  1*. 

V(eia)  -  VQ.ix  )  ♦  0.25  [VC,  )  +  V( ,,4  )].  (16) 

From  Eqs.  (3d): 

V(*»a  )  -  V(,l,  ),  (17) 

V(*a  «)-▼(••  4>-  (18) 

THE  STRESS  TENSOR 


We  define  Q  ■  S*/(3Sx  +•  1/2S,  )  (which  le  [”v/l-2v]  for  »n.  isotropic 
solid).  Then  Eq.  (4)  may  be  written,  for  the  diagonal  stress  components,  as: 

j  -  (1/2S1)_:L[(1-Q)ell  -  Qekk  -  Qe^].  (19) 

Here  1  -  1,2,3;  j  -  2,3, l;k  -  3,1,2. 

From  Eq.  (19): 

V(0ii)%  "  d'2**)'1  ((l-Q)2V(eii)  +  Q2[V(ekk)  +  V(«j j )]}^.  (20) 

Therefore,  with  Eqs.  (13-15): 

V(<^x  )^  -  (1/2S,  )-1  {(2-4Q  +  4 C?)V(,fx)  +  <?V(t,  fx  ) 

+  (1-2Q  +  3<f  V(I)}*,  (21) 

V(q,a  )*  -  (l/2Sj  )_1  {(1-2Q  +  4<?  )V(,  ^  )  +  (1-2Q  +  <?  )V(,,  £x  ) 

♦  (1-2Q  +  3<?)V(I)}*,  (22) 

V(q,,)%  -  (1/2S,  )"1{(1-2Q  +  4qf  )V(,  lx  )  ♦  <fv(,,lx) 

♦  (1-2Q  +  3<f  )V(I)}% 


(23) 


Similarly: 


V(oia)%  -  (1/2%  )"1[V(4S  %  )  +  0.25  [V( ,%)  +  V(t0%)]*,  (24) 

V(ql,  )%  -  (1/2%  )_1V(o  %  )*,  (25) 

V(%,  )  -  (1-2%)‘1V(,0%)%.  (26) 

EXAMPLES 

To  illustrate  the  typical  magnitudes  of  the  errors  due  to  counting 
statistics,  we  employed  data  from  Ref.  5,  for  a  ground  steel  specimen,  that 
is  we  used  the  peak  positions  and  the  variances  in  these  positions  with  Eq. 
(9).  [Formulae  to  calculate  this  variance  for  the  parabolic  fit  employed 
in  Ref.  5  are  given  in  Ref.  1;  for  other  types  of  fits  the  appropriate  eq¬ 
uation  may  be  substituted.]  The  resultant  errors  are  given  in  Tables  I-III. 
For  the  first  two  tables  it  was  assumed  that  the  error  in  d.  .  was  the 
actual  measured  value.  If  there  is  no  preferred  orientation,  the  intensity 
of  the  peak  changes  little  with  the  ^-tilt.  In  this  case.  Tables  I  and  II 
show  the  effect  of  the  uncertainty  in  % ;  reducing  this  error  all  the 
stress  components  by  the  same  proportion,  except  o^,  ,  ,  which  remain 

relatively  unaffected,  because  thd  role  of  the  error  in  do  is  damped  by 
(d-P  in  Eq.  (12). 

If  there  is  preferred  orientation,  the  peak  intensity  can  vary  greatly 
with  i>  and  there  will  be  large  variances  contributing  to  V(%  )  from  weak 
peaks.  This  was  minimized  in  the  following  way.  The  average  variance,  o^, 
in  the  26  peak  position  for,-*^  and  was  obtained  and  the  weighting  factor 
c^  was  formed: 

Ci  -  (l/o^)  /  l  (l/o^)  (27) 

The  Eqs.  11  and  12  were  then  altered  to  multiply  Vi(26+),  ^(26*) 
terms  by  this  weighting  for  Table  III.  There  is  only  a  small  difference 
(between  Tables  II  and  III)  because  of  the  lack  of  texture  in  the  specimen; 
the  peak  intensity  changed  only  by  about  8  pet  with  <p.  With  more  severe 
preferred  orientation  the  effect  will  be  larger. 

CORCLUDIftG  REMAEES 


There  are  now  adequate  equations  for  calculating  errors  in  stress 


TABLE  IX:  STRESS  TENSOR  AND  STANDARD  DEVIATIONS  WHEN 


V(d,  >*  -  0 

.00004  A  * 

539.74 

(48.24) 

DATA  SET  1 
-24.03  (24.50) 

552.16  (47.26) 

-39.15 

2.30 

80.41 

(4.58) 

(3.56) 

(38.96) 

520.60 

(40.52) 

DATA  SET  2 
-4.03  (19.95) 

555.19  (39.73) 

-34.17 

0.11 

82.20 

(3.21) 

(2.69) 

(32.75) 

535.03 

(47.81) 

DATA  SET  3 
-20.14  (24.35) 

555.98  (46.81) 

-40.19 

-0.98 

36.66 

(5.72) 

(4.56) 

(38.84) 

538.53 

(42.84) 

DATA  SET  4 
-30.63  (21.10) 

565.37  (42.51) 

-38.03 

0.76 

88.18 

(3.83) 

(3.89) 

(34.67) 

♦values 

given  In  MPa 

TABLE  III:  WEIGHTED 

STRESS  TENSOR  AND  STANDARD  DEVIATIONS* 

536.24 

(48.56) 

DATA  SET  1 
-24.62  (24.56) 

554.65  (47.60) 

-38.33 

2.90 

80.68 

(4.59) 

(3.65) 

(39.22) 

520.29 

(41.84) 

DATA  SET  2 
-6.04  (20.11) 

560.22  (40.60) 

-34.82 

1.56 

85.29 

(3.55) 

(2.73) 

(33.76) 

532.56 

(48.15) 

DATA  SET  3 
-7.90  (24.93) 

549.83  (47.16) 

-39.66 

-2.81 

82.21 

(5.81) 

(4.57) 

(39.16) 

539.23 

(42.88) 

DATA  SET  4 
-31.22  (21.23) 

565.17  (42.65) 

-38.28 

-0.59 

88.98 

(3.85) 

(3.94) 

(34.68) 

* 

< 

* 

*  -  0.00004  A; 

values  given  in  MPa 

measurements  due  to  instrumental  effects,  counting  statistics  and  in  the 
x-ray  elastic  constants.  We  would  like  to  conclude  this  paper  with  a  plea 
to  the  community  making  stress  measurements  via  diffraction  to  regularly 
report  these  errors  with  their  data.  It  is  all  too  common  for  investigators 
to  repeat  a  measurement  (of  stress  or  an  elastic  constant)  once  and  to  use 
the  difference  as  an  error  estimate.  Another  practice  is  to  dust  a  stress- 
free  powder  on  the  specimen  surface  and  to  use  a  (single)  measurement  of 
the  stresses  measured  with  this  powder  as  an  error  estimate.  Finally,  some 
report  an  error  in  a  slope  vs.  sit^4>  obtained  by  least- squares,  but  ignore 
the  uncertainty  in  each  point  in  this  plot  in  estimating  errors.  None  of 
these  procedures  is  particularly  satisfying  in  a  statistical  sense.  Of 
course,  if  time  permits,  the  average  of,  say,  ten  repetitions  of  a  measur¬ 
ement  is  the  best  of  all  error  estimates.  If  this  cannot  be  done,  error 
estimates  from  the  available  equations  are  far  more  satisfactory  than  the 
currently  all  -  too  common  procedures. 
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